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Abstract 

We show how the superintegrability of certain systems can be deduced from the presence of 
multiple parameters in the rational Lax matrix representation. This is also related to the fact 
that such systems admit a separation of variables in parametric families of coordinate systems. 



1 Rational Lax matrix representations of integrable 

systems 

1.1 Classical i?-matrix theory of commuting isospectral flows 

In the classical i?-matrix approach to finite dimensional integrable systems |7J 121 E] , there 
is a Poisson map from the phase space into a space of r x r Lax matrices M{\) depending 
rationally, trigonometrically or elliptically on a spectral parameter A. The Poisson bracket 
is denned by the relation 

{yV(A) fJVO)} = [r(A - N{\) ® I + 1 ® N{n)) , (1) 



1 Work supported in part by the Natural Sciences and Engineering Research Council of Canada 
(NSERC) and the Fonds FCAR du Quebec. 
2 harnad@crm . umontreal. ca 
3 yermolae@crm.umontreal.ca 



1 



where both sides are interpreted as elements of End(C r ®C r ). The symbol { ® } signifies a 
simultaneous tensor product in End(C r ) ~ gl(r) and Poisson brackets in the components, 
and r(A — fi) denotes the classical R— matrix. The simplest case is the rational i?-matrix, 




P i2 (u <g> v) := v ® u , 



(2) 



with jV(A) a rational function of A. 



A/-(A) 



(3) 



i=l a=l 



(A - 



B(X) 



7V ia , 5i G gl(r) . 



(4) 



Equations (Q), (j2J) define the standard linear, rational R- matrix structure. 

It follows from the properties of classical R- matrices j3 E] that elements of the algebra 
of spectral invariants <fi(N) G X(gl(r)) Poisson commute amongst themselves and generate 
commuting isospectral flows determined by the Lax equations: 



where M is here thought of as an element of the loop algebra gl(r), identified in a standard 
way with its dual gl (r) through the trace-residue pairing, and ( . )± denotes projection 
to the ± components relative to the usual splitting of the loop algebra into positive and 
negative components 



(i.e. those admitting holomorphic continuations to the interior (+) and exterior (-) of the 
unit circle respectively with the latter normalized to vanish at oo). The spectral invariants 
generate a maximal Poisson commuting algebra on generic symplectic leaves, defining 
completely integrable systems E]; i.e., there are as many functionally independent 
generators as half the dimension of the leaf. 

1.2 2x2 rational Lax matrices 

In the following, we shall limit our discussion to the case of 2 x 2 Lax matrices, although 
most of the considerations that follow are easily extended to higher rank. We may without 
loss of generality take M{\) to be traceless (since the trace coefficients are Casimirs) 



Af=±[(dM ± ,Af] , 



(5) 



gl(r)) = gl(r))+ + gl(r)) 



(6) 




(7) 
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where the rational functions e(A), /(A), h(X) satisfy the Poisson bracket relations 



(MA), e(„)} = 

(MA), M) = -^4^ 

{ e(A), /(„)} = _ 2 M^M . (8 ) 

For this case, the ringX(gl(2)) of spectral invariants, when restricted to the symplectic 
leaves of the i?-matrix Poisson structure, is generated by the quadratic trace invariants; 
i.e., the coefficients determining the numerator of the rational function 

A(A) := -itr(A/" 2 (A)) = h\\) - \ (e(A)/(A) + /(A)e(A)) . (9) 

(The order in the last two terms is irrelevant of course, but it is written here in a form 
that will also be valid in the quantum version below.) If, for example, the polynomial 
part £>(A) of jV(A) is taken to vanish, and only first order poles appear in jV(A), we have 

i=i % 

i=l % 

W> - do) 
i=l A a% 

where the quantities {e^, fi,hi}i 

=l...n 8X6 db set of n sl(2) generators, which may be canon- 

ically coordinatized as: 

f _ 1 2 

hi := -Xiiji, i = l,...n, (11) 

where {/x?}i=i... n are the values of the sl(2) Casimir invariants and {xi, yi}i=i... n form a 
set of canonical coordinates on the symplectic leaves . 
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1.2.1 Parametric dependence of invariants and superintegrability 

Again, taking the case when the polynomial part B(\) of jV(A) vanishes (but not neces- 
sarily just first order poles), a complete set of generators is given by 

(j) ta :=res A = Ql (A-a i ) a tr(Ar 2 (A)), i = 1, . . .n, a = 0, . . .n* - 1 . (12) 

These commute amongst themselves, but they each depend upon the pole locations 
{cti}i=i... n in jV(A). However, the linear combination: 

n 

4>si := ^2 ai(fiio = res A =ootr(A/' 2 (A)) (13) 

i=l 

does not depend on the c^'s. In general, there is no reason for the invariants 0i O («i) to 
commute with each other for different choices of the a^s. But, regardless of the values 
chosen, they will commute with (f>si- Since the </>j a (a!j)'s for different choices of a^s in 
general do not generate the same algebra of functions, we may conclude that, taken 
together, for different evaluations of the parameters there are more functionally 

independent integrals that Poisson commute with <f>si than half the dimension of the 
symplectic leaf, and hence the Hamiltonian system it generates is superintegrable. (In 
fact, in most cases, it may be shown to be maximally superintegrable; see the examples 
below.) 

In particular, if we take the case of purely simple poles as above in (jlOj). the resulting 
Hamiltonian is: 

i=l j=l i=l i=l i=l 1 

which, when constrained to the (co)tangent bundle of the n — 1 sphere S n ~ 1 

n n 

X^ 2 = 1 ' 53^ = 0, (15) 

i=l i=l 

yields the superintegrable system 

"- = iE^ + iE§. < 16 > 

j=l i=l * 

which is the trivial case of the Rosochatius system (without a harmonic oscillator poten- 
tial). 
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1.2.2 Separation of variables 



Another viewpoint that helps to explain the superintegrability of systems arising in this 
way is to note that they may be completely separated in a canonical coordinate system 
determined by the values of the pole parameters which, for the sl(2) case with simple 
poles with the phase space constrained to S 1 ™ -1 as above reduces to the sphero-conical 
system {A*, (i}i=i...n-i defined by: 



, A • <>; nr=i(A"« 4 )' 2^{X-a i ) ' 1 ) 

These are just the points (Aj, Q) on the invariant spectral curve 

C 2 + ^A(A) = (18) 

where the matrix element /(A) vanishes and Q = h(\) are the eigenvalues at these points. 
These are particular cases of the spectral Darboux coordinates of [3J E] . (Note that these 
become hyperellipsoidal coordinates if there is a constant term added in the definition 

flUB of /(a).) 

The point to note is that the separation of variables occurs in these coordinates si- 
multaneously for all the invariants <fii a , viewed as generators of Hamiltonian flows. But 
again, since the leading term spectral invariant <ftsi does not depend on the values of the 
parameters c^, it admits a separation of variables in any of the family of sphero-conical 
(or hyperellipsoidal) coordinates obtained by varying these parameters. This simultane- 
ous separability in multiple coordinates may be viewed as an alternative explanation of 
the origin of the superintegrability of such systems. (In fact, both these viewpoints are a 
result of the classical r- matrix setting, and in a sense may be considered as equivalent.) 

In the examples given below in the following section, the same principle is used to 
deduce superintegrable systems from sl(2) Lax matrices satisfying the Poisson bracket 
relations (IU). 



1.2.3 Quantum integrable systems 

The above discussion is easily extended to the canonically quantized version of such sys- 
tems. All that must be done is to replace the matrix elements defining J\f(X) by their 
quantized forms e(A), /(A), h(X), which must satisfy the commutator analogs of the Pois- 
son bracket relations (jHJ) 

[h{\), e{ji)\ = r 
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(MX), /</.)] = - f MzM 

A — /Ji 

m, /ooi = -M=M, . (19) 

A — 11 

These can be realized by canonical quantization of the underlying classical phase space 
variables. For example, in the case of simple poles only, with vanishing polynomial term 
B(X), we have: 



e(A) := £ 



A - a* 

n 



/(A) := £ /; 



A — OLi 

i=l 

hi 



h(X) := £y^, ( 2 °) 
where the sl(2) generators {e^, /j, ^} may be represented by the operators 

d 2 nT 



2 \dxj x\ 

' = 2, Xi 

1/9 1 



:= 2^^ + 2j' < = 1 .- n . ( 21 ) 

and the commuting invariants are similarly given by the coefficients of the numerator 
polynomial of the quantum spectral invariant: 



A(A) := h\\) - \ (e(A)/(A) + /(A)e(A)) . (22) 



The resulting systems are similarly quantum integrable, and separable in the same 
coordinates as the classical ones jH] and, for the same reasons as above, the quantum 
version of the Hamiltonian (f>si is superintegrable. 

In the following section, a number of examples of such classical and quantum super- 
integrable systems will be given. 



2 Examples of superintegrable classical and quantum 
systems 

The examples given below arise in the framework of the so-called Krall-Scheffer problem 
jH] of describing all two-dimensional analogs of classical orthogonal polynomials which 
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result in nine classes of second-order partial differential equations on the plane or on 
constant curvature surfaces. It was shown in that all nine cases are connected with 

superintegrable systems. The following are some illustrative examples. 



2.1 Example 1. The sphere. 

2.1.1 Classical Lax Matrix 

The first case corresponds to three simple poles and vanishing B(X) . The Lax matrix 
has the form: 

A rm_ ^ i N * i _(h{\) /(A) \ 

^ (A) - (A^y + (X^y + (X^y - 1 e(A) -ma) J (23) 

where the matrix elements of the iVj generate a Poisson bracket realization of (sl(2)) 3 : 

N if «, p. 2 + % , (24) 

iV 2 = - f S2P l + $ ) (25) 

N 3 = i ( s « *»' + f ) (26) 

Here (pi,P2,Ps) are canonically conjugate to (si,s 2 ,s 3 ) (and these coincide with the co- 
ordinates {xi, Ui}i=i... n above). 

2.1.2 Commuting invariants 

The invariants are the coefficients of: 

A + ^ + A + ^ + + (27) 

Note that only two of the integrals Hi, H 2 and if 3 are independent, since their sum 
is zero. The Hamiltonian of the problem is given by their linear combination: 



H = aH 1+ (3H 2 + 7 tf 3 = k P i 2 + P2 2 + p 3 2 ) + ^ + ~2 + ~2~ & 

I S\ S2 S3 
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This describes the Rosochatius system with harmonic oscillator terms absent on the cotan- 
gent bundle of a two- sphere in M 3 : 

si 2 + s 2 + s 3 2 = 1, sipi + s 2 p 2 + s 3 p 3 = . (29) 

The integrals Hi, H 2 and if 3 follows: 



Hi 



1 L 13 2 + /i 3 2 si 2 /s 3 2 + Vi 2 s 3 2 /si 2 1 L 12 2 + fi 1 2 s 2 2 /s 1 2 + fi 2 2 s 1 2 /s 2 2 



2 ot-7 2 a-/3 

1 L 23 2 + fX 3 2 S 2 2 /s 3 2 + /^V/V , 1^12 2 +^l 2 S2 2 /Sl 2 + /i 2 2 S 1 2 /s 2 2 



2 2 0-7 2 a-/3 

1 L 23 2 + /i 3 2 S 2 2 /s 3 2 + H2 2 S 3 2 /S 2 2 , 1 L 13 2 + fi 3 2 Sl 2 /s 3 2 + m 2 S 3 2 /.Si 2 

2 p — 7 z a — 7 

where = «ip 2 — S2P1 . 

Note that the Hamiltonian H is independent of the parameters (a,/3, 7), whereas the 
invariants iii, if 2 do depend on them. Therefore, different choices for the parameters give 
distinct integrals that commute with H, but do not commute with each other. 



2.1.3 Separating coordinates 

The separating coordinates (Ai,A 2 ) in this case are sphero-conical coordinates. The cor- 
responding momenta are denoted (£i,£ 2 ) • They are related to (si,s 2 ,s 3 ) and (pi,p 2 ,p 3 ) 
by: 

2 _ (a- Ai)(a- A 2 ) I, S1P1 s 2 p 2 . -siPi - s 2 p 2 

{a — p){a — r y) 2 Ai — a Ai — p Ai — 7 

2 (/?- Ai)(/3- A 2 ) 1, sipi s 2 p 2 -sipi-s 2 p 2 

(/?-a)(/?-7) 2 A 2 - a X 2 - p A 2 - 7 

(7 - a) (7-/?) 

2.1.4 Quantum system 

The quantum versions of the integrals above, denoted Hi, H 2 , H 3 , are obtained by replac- 
ing the matrix elements of N(X) by the corresponding differential operators, e(A), /(A), h(X), 
which in the case of simple poles are as in (J2(J|) - (J22|) . 

The quantization procedure leads to replacing the Lj/s by their quantum version: 



Lij = V-l(sid/dsj - sjd/dsi) (34) 
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Introducing the functions 

^jk '■= Vj 2 Sk 2 /sj 2 + l^k 2 Sj 2 /s k 2 j, k = 1..3 (35) 
and denoting a = ai , f3 = a 2 , 7 = 0:3, we can present the quantum integrals as 

l£^±i4 a = 1..3 (36) 



The quantum Hamiltonian is 



H = -hd, 2 + d 2 2 + d 3 2 ) + ^ + 4 + — 2 - (37) 

I S\ S2 S3 



The separating coordinates are the configuration space part of the ones for the classical 
case (Ai, A 2 ). 

2.2 Example 2. The hyperboloid. 
2.2.1 Classical Lax Matrix 

Consider now a Lax matrix with one first order and one second order pole: 

iVi No N 3 

N ^ = J^ + (x^ + (^Wy < 38 > 

with 

1 ( + S2P2 2pip 2 + 27x72 



2 \ -2s!S 2 -S1P1 - S2V2 

N = \ -v, 2 -il 

2 V S 2 S 2 P! 

N 3 = l -( S * P \ P^ + ^l) (39) 
2 V S3 2 -S3P3 J y 1 

Here we have introduced the following notations 

27172 := — 3 -2— , 7 2 := --r , 7 3 := -3 ■ (40) 

5 2 o 2 * 2 S3 

The matrix elements of (Ni,N 2 ) generate a Poisson bracket realization of the jet 
extension sl(2)^^* while those of N 3 generate a second sl(2). 
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2.2.2 Commuting invariants 

The trace formula again gives us only two independent commuting invariants Hi and H 2 

1 N (\\2 = -^1 #2 ^1/^2 /^2 2 #3 f^3 2 / 41 x 

2 lJ (A -a) (A -a) 2 (A - a) 3 2(A - a) 4 (A - /3) 2(A-/3) 21 j 
since, by taking the residue we obtain 

Hi + H 3 = . (42) 
The superintegrable Hamiltonian in this case is: 

H = (a- (3)Hi + H 2 - i/i 3 2 = 2 Pl p 2 - Pl 2 + p 3 2 + 2 7l72 - 7s 2 + 7 f . (43) 
The quadratic constraint now defines a hyperboloid 

2sis 2 + s 3 2 = 1 • (44) 
In the ambient coordinates the integrals Hi and H 2 are 

H (sm - s 3 p 2 )(s 3 pi ~ s 2 p 3 ) - 7 fgis 2 - 2 7l72 s 2 . _ ((s 3 pi - s 2 p 3 ) 2 + 7 fs 2 + 7 2 s 2 

1 a -(3 2{a-(3) 2 

#2 = -(S1P1-S2P2) -2 7 i 72 si + 2 _ ^ . (45) 

Again, whereas the Hamiltonian H does not depend on the parameters (a, (3) the 
integrals Hi, H 2 do, which provides an explanation for the superintegrability in this case. 

2.2.3 Separating coordinates 

These are determined by the relations: 

2 _ (ff- A 2 ) p \, sipi + S2P2 S2P1 ^3P3 v 

53 ~ (a-/5) 2 ^ 2 l A x -a (A x - a) 2 + A x - /3 ] 1 j 



2 = (a - Ai)(a - A 2 ) + S2P2 S2P1 s 3 p 3 

32 (a-(3) ^ 2 2 l A 2 -a (A 2 - a) 2 + A 2 - (3 } 1 } 

SlS2 ""2 ( (a~^ 1} - (48) 
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2.2.4 Quantum system 

The quantized integrals Hi,H 2 , H 3 are obtained as before by replacing all conjugate vari- 
ables by the corresponding differential operators. The quantum integrals may then be 
expressed as 

rj Qid 3 - s 3 <9 2 ) (s 3 di - s 2 d 3 ) + 7fsis 2 + 27172s 2 , (s 3 di - s 2 <9 3 ) 2 - ^s 2 , - 7|s§ 

1 a-P 2(a-(3) 2 

H 2 = -L l2 -2 lll2Sl ^-—^ , (49) 

where dk :— d/dsk- 

The quantum Hamiltonian is 

H = 2<9k9 2 - d± 2 + d 3 2 + 2 7l72 - 7 2 + 73 2 , (50) 

and this again separates in the configuration space coordinates (Ai, A 2 ). 

2.3 Example 3. The plane. 
2.3.1 Classical Lax Matrix 

For the cases with zero curvature like the example to follow, the polynomial part £>(A) 
of the Lax matrix does not vanish. The simplest case involves two distinct finite poles in 
iV(A) and constant B(X) 



N(X) 



—a \ 1 / sipi pi 2 + - 1 



2 
1 

72" 



10 ; 2(A-a) ^ - Sl 2 - SlPl 
i 1 ( s 2 p 2 p 2 2 + 4\ ( U 

+ ^W){s 2 i - S2P ; ) • (51) 

The matrix elements of the residues N ± , N 2 generate two copies of sl(2). 

2.3.2 Commuting invariants 

The invariants of motion are defined by: 
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The superintegrable Hamiltonian in this case is given by 

H = ^re Soo triV(A) 2 = ]( Pl 2 + p 2 2 + a( Sl 2 + s 2 2 ) + ^ + ^! 
4 4 Si z S2 



which gives an isotropic oscillator together with Rosochatius terms. As before (pi,p 2 ) are 
canonically conjugate to (si, s 2 ). 

In terms of the ambient space coordinates the integrals H\ and H 2 are : 



tt 2 , „ n 2 , A*l ' / 7" 2 I A t l S 2 | / i 2 S l\ 

^ = Pi+asi+ 7?-2(^) (Li2+ ^r + ^r ) 



A*2 , 1 / r 2 | /flfl i A*2 S 1' 

S | + 2(a-/3) (Ll2+ 



tf 2 = ^ + a S 2 -^ + — ^— (L 2 2 + ^ + ^i). (54) 



Where L 12 := S1P2 _ S2P1 and = + H2). Here the additional integral results from 
the parametric dependence on (a — (3). 

2.3.3 Separating coordinates 

The separating coordinates (Ai, A2,^i,^) in this case are defined by 

51 = 2 — (^ — 6 = "2 ( a7^ + a7^ } (55) 

2 (a- Ai)(a- A 2 ) 1, S2P2 v , R v 

52 ~~ 2 (^9) 6 -"2 ( A 2 ^ + A^ ) (56) 

2.3.4 Quantum system 

The Hamiltonian of the corresponding quantum problem is 

H = ire Soo triV(A) 2 = -(d^ + d 2 2 + a( Sl 2 + s 2 2 ) + ^ + ^4) = 7(^1 + #2). (57) 
4 4 S2 4 

The quantum integrals ifi and H 2 are: 

^ - 9i+asi+ ^~2(a^) (Li2+ ^r + ^r ) 

2 1 ,,2„2 ,,2 2 



= % + as *-% + ^f L » + t f + "f ' <58) 

and the separating coordinates are again (Ai, A 2 ), which depend on the additional param- 
eter (a — (3). 



12 



References 



[1] M.R. Adams, J. Harnad and J. Hurtubise, "Darboux Coordinates and Liouville- 
Arnold Integration in Loop Algebras", Commun. Math. Phys. 155, 385-413 (1993). 

[2] L.D. Faddeev and L.A. Takhtajan, Hamiltonian Methods in the Theory of Solitons, 
Springer- Verlag, Berlin (1986). 

[3] J. Harnad, "Isospectral Flow and Liouville-Arnold Integration in Loop Algebras=D3, 
in: Geometric and Quantum Methods in Integrable systems, Springer Lecture Notes 
in Physics 424, ed. G. F. Helminck (Springer- Verlag, N.Y., Heidelberg, 1993). 

[4] J. Harnad, "Loop Groups, R-Matrices and Separation of Variables" , In: Integrable 
Systems: From Classical to Quantum , CRM Proceedings and Lecture Notes Series 
26, 21- 54. (AMS, Providence, R.I, 2000). 

[5] J. Harnad, L. Vinet, O. Yermolaeva and A. Zhedanov, "Two-dimensional Krall- 
Sheffer polynomials and integrable systems", J. Phys. A : Math. Gen. 34, 1-7 (2001). 

[6] J. Harnad and P. Winternitz, "Classical and Quantum Integrable Systems in Lgl(2) + 
and Separation of Variables" , Commun. Math. Phys. 192, 263-285 (1995). 

[7] Semenov-Tian-Shansky, M. A., "What is a Classical R-Matrix", Fund. Anal. Appl. 
17, 259-272 (1983) 

[8] H. Krall, I. M. Sheffer, "Orthogonal polynomials in two variables", Ann. Math. Pure 
Appl.,, V.76, 325-376 (1967). 

[9] J. Harnad, O. Yermolaeva and A. Zhedanov, "R-matrix approach to the Krall-Sheffer 
problem", Proc. Inst. Math. 43, Part 1, 314-321, (2001), (NAS, Ukraine), eds. A.G. 
Nikitin, V.M. Boyko and R.O. Popov. 



13 



